Precision Measurements of the Nucleon Strange Form Factors at ~ 0.1 GeV^ 
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We report new measurements of the parity-violating asymmetry ^pv in elastic scattering of 
3 GeV electrons off hydrogen and ^Ife targets with {Oiat) ~ 6.0°. The ''He result is Apv = (+6.40 ± 
0.23 (stat) ±0.12 (syst)) x 10"^. The hydrogen resuh is Apv = (-1.58±0.12 (stat) ±0.04 (syst)) x 
10~^. These results significantly improve constraints on the electric and magnetic strange form 
factors Ge and Gm. We extract G% = 0.002 ± 0.014 ± 0.007 at {Q^) = 0.077 GeV^ and Ge + 
0.09 Gm = 0.007 ± 0.011 ± 0.006 at (Q^) = 0.109 GeV^, providing new limits on the role of strange 
quarks in the nucleon charge and magnetization distributions. 

PACS numbers: 25.30.Bf, 13.60.Fz, ll.30.Er, 13.40.Gp, 14.20.Dh 



Over the past several decades, high-energy lepton- 
nucleon scattering has revealed the rich structure of 
the nucleon over a wide range of length scales. In re- 
cent years, increasingly sensitive measurements of elas- 



tic electron-nucleon scattering, mediated by photon ex- 
change and Z" exchange, have enabled the measurement 
of the electromagnetic and neutral weak form factors. 
These functions of the 4-momentum transfer charac- 
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terize nuclcon charge and magnetization distributions. 

In particular, the neutral weak form factor measure- 
ments provide a way to probe dynamics of the "sea" 
of virtual light (up, down and strange) quark-antiquark 
pairs that surrounds each valence quark in the nucleon. 
Since the TP boson couples to various quarks with differ- 
ent relative strengths compared to the photon, a com- 
bined analysis of proton and neutron electromagnetic 
form factor and proton neutral weak form factor mea- 
surements, along with the assumption of charge symme- 
try, allows the determination of the strange electric and 
magnetic form factors G% and G\j [1, 2]. 

The established experimental technique to mea- 
sure the electron-nucleon weak neutral current am- 
plitude is parity- violating electron scattering [3, 4]. 
Longitudinally-polarized electron-scattering off unpolar- 
ized targets can access a parity-violating asymmetry 
Apv = {a-R- (Tl)/{(tr + (Tl), where crfl(i) is the cross 
section for incident right(lcft)-handcd electrons. Arising 
from the interference of the weak and electromagnetic 
amplitudes, Apy increases with [5]. 

Four experimental programs have been designed to ac- 
cess the range of 0.1 to 1 GeV^, where the Apy expec- 
tations range from one to tens of parts per million (ppm). 
The published measurements [6 12] arc mutually consis- 
tent. An intriguing pattern in the low-Q^ behavior seen 
in [9, 10] has marginal statistical significance. 

In this paper, we significantly improve our two previous 
measurements [11, 12] of Apy in elastic electron scatter- 
ing from and ^He nuclei. Since ^pv for is sensitive 
to a linear combination of G|; and while that for ''He 
is sensitive only to G|;, a simultaneous analysis of both 
measurements results in the most precise determination 
to date of G% and Glj at _ 0.1 GcV^. 

The measurements were carried out in Hall A at the 
Thomas Jefferson National Accelerator Facility (JLab). 
As described in detail in two previous publications [11, 
12], a 35 to 55 fiA continuous-wave beam of ~3 GeV lon- 
gitudinally polarized electrons was incident on 20 cm long 
cryogenic targets. Elastically scattered electrons were fo- 
cused into background-free regions by a symmetric pair of 
high-resolution spectrometer systems. The scattered flux 
was intercepted by identical detector segments in each 
arm (two for ^H, one for "^He), resulting in Cherenkov 
light collected by photomultiplier tubes (PMTs). 

The hclicity of the electron beam, generated by pho- 
toemission off a GaAs wafer, is determined by the hand- 
edness of the incident laser light's circular polarization. 
This was selected pscudorandomly at 15 Hz and toggled 
to the opposing helicity after 33.3 ms, with each of these 
equal periods of constant helicity referred to as a "win- 
dow." PMT and beam monitor responses for two consec- 
utive windows of opposite helicity were integrated, digi- 
tized, and grouped as a "pair" for asymmetry analysis. 

The beam monitors, target, detector components, elec- 
tronics and accelerator tune were optimized such that the 





Helium 


Hydrogen 


-^Intensity 


-0.377 ppm 


0.406 ppm 


■^Energy 


3 ppb 


0.2 ppb 


Ax 


-0.2 nm 


0.5 nm 


Ax' 


4.4 nrad 


-0.2 nrad 


Ay 


-26 nm 


1.7 nm 


Ay' 


-4.4 nrad 


0.2 nrad 



TABLE I: Average beam asymmetries under polarization re- 
versal in intensity and energy and differences in horizontal 
and vertical position (Aa;, Ay) and angle (Aa;', Ay') . 

fluctuation in the PMT response over a pair was domi- 
nated by counting statistics of the scattered flux for rates 
up to 100 MHz. This facilitated ^pv measurements with 
statistical uncertainty as small as 100 parts per billion 
(ppb) in a reasonable length of time. To keep spurious 
beam-induced asymmetries under control at this level, 
the laser optics leading to the photocathodc were care- 
fully designed and monitored. Indeed, averaged over the 
entire period of data collection with the hydrogen tar- 
get, the achieved level of control surpassed all previous 
benchmarks, as summarized in Table I. 

The data collection took place over 55 days ("^He) and 
36 days (^H). A half-wave (A/2) plate was periodically 
inserted into the laser optical path which passively re- 
versed the sign of the electron beam polarization. With 
roughly equal statistics in each state, many systematic 
effects were suppressed. There were 121 (^He) and 41 
(^H) such reversals. The data set between two successive 
A/2 reversals is referred to as a "slug." 

Loose requirements were imposed on beam quality to 
remove periods of instability, leaving about 95% of the 
data sample for further analysis. No helicity-dcpendcnt 
cuts were applied. The final data sample consisted of 
35.0 X 10*^ (*He) and 26.4 x 10*^ (^H) pairs. The right-left 
helicity asymmetry in the integrated detector response, 
normalized to the beam intensity, was computed for each 
pair to form the raw asymmetry A^aw ■ The dependence of 
Aj-aw on fluctuations in the five correlated beam param- 
eter differences Axi is quantified as ^beam = X^c^Axi, 
where the coefficients Ci quantify the ^raw beam parame- 
ter sensitivity. The clcctroweak physics of the signal and 
backgrounds is contained in A^orr = ^raw — ^beam- 

The AcoTT window-pair distributions for the two com- 
plete data samples were perfectly Gaussian over more 
than 4 orders of magnitude with RMS widths of 1130 
ppm (*He) and 540 ppm (^H); the dominant source of 
noise in the PMT response was counting statistics. To 
further test that the data behaved statistically and the er- 
rors were being accurately calculated, Acott averages and 
statistical errors for typical one hour runs, consisting of 
about 50k pairs each, were studied. Each set of roughly 
400 average ^corr values, normalized by the correspond- 
ing statistical errors, populated a Gaussian distribution 
of unit variance as expected. 
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(DOF = 1 
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A 


4.80±0.27 
5.12±0.27 


0.75 
0.78 


-5.41±0.27 
-5.38±0.27 


1.12 
1.07 


5.10±0.19 
5.25±0.19 


0.95 
0.92 




(DOF = 20) 


(DOF = 19) 


(DOF = 40) 


^corr 


-1.40±0.15 
-1.41±0.15 


0.73 
0.81 


1.42±0.15 
1.43±0.15 


1.04 
1.02 


-1.41±0.11 
-1.42±0.11 


0.86 
0.89 



TABLE II: Raw and corrected asymmetries (in ppm) and re- 
duced "slug" (rx^)) broken up by A/2 reversals. The dif- 
ferences between Aram and Acott result from corrections for 
energy, position, and angle differences which are summarized 
in Table I. 



Systematic effects in Abeam estimations were studied. 
When averaged over all detector segments, the coeffi- 
cients Ci were much smaller than those for individual de- 
tector segments due to the symmetric geometry of the 
apparatus. Limits on systematic uncertainties in the Ci's 
in the range of 10 to 30% were set by inspecting resid- 
ual correlations of ylcorr's of individual detector segments 
with helicity-correlated beam asymmetries. 

Another important validation was to use two indepen- 
dent methods to calculate q. The first relied on linear 
regression of the observed response of the detector PMTs 
to intrinsic beam fluctuations. The other used calibration 
data in which the beam was modulated, by amounts large 
compared to intrinsic beam fluctuations, using steering 
magnets and an accelerating cavity. Differences in the 
two Abeam Calculations were always much smaller than 
corresponding Acorr statistical errors. 

Final Aeon results were calculated using the beam 
modulation technique and are summarized in Table IL 
Due to the excellent control of beam parameter differ- 
ences Axi summarized in Table I, Aeon- — Araw values are 
of the order of, or much smaller than, the corresponding 
statistical errors. Under A/2 reversal, the absolute values 
of Aeon are consistent within statistical errors. The re- 
duced for Aeon "slug" averages is close to one in every 
case, indicating that any residual beam-related system- 
atic effects were small and randomized over the time pe- 
riod of A/2 reversals (typically 5 to 10 hours). The final 
Aeon results are A^„°, = +5.25 ± 0.19(stat) ± 0.05(syst) 
ppm and A^,, ^ -1.42 ± O.ll(stat) ± 0.02(syst) ppm. 

The physics asymmetry Aphys is formed from Aeon, 

. _ K Aeon - Ph J2i ^ifi /in 

with corrections for the beam polarization Pb, back- 
ground fractions fi with asymmetries Ai and finite kine- 
matic acceptance K. These corrections are described be- 
low and summarized in Table III. The first line lists the 
cumulative Abeam corrections discussed above, scaled by 
if/Pb. 

A powerful feature of the apparatus is the spectrom- 
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FIG. 1: Single-particle spectra obtained in dedicated low- 
current runs. The insets show the same spectra on a loga- 
rithmic scale. The vertical lines delineate the extent of the 
detectors. Inelastic scattering from *IIe is entirely contained 
in the hatched area. The shaded regions, visible only in the 
log plots, show the contribution from target windows. 



Correction (ppb) 


Helium 


Hydrogen 


Beam Asyms. 
Target window bkg. 
Helium QE bkg. 
Rescatter bkg. 
Nonlinearity 


183 ± 59 
113 ± 32 
12 ± 20 
20 ± 15 
± 58 


-10 ± 17 
7 ± 19 

2 ± 4 
± 15 


Scale Factor 


Helium 


Hydrogen 


Acceptance factor K 
Q2 Scale 
Polarization Pb 


1.000 ± 0.001 
1.000 ± 0.009 
0.844 ± 0.008 


0.979 ± 0.002 
1.000 ± 0.017 
0.871 ± 0.009 



TABLE HI: Corrections to Aeon and systematic errors. 



eters' ability to focus the elastically scattered electrons 
into a compact region. Indeed, much less than 1% of 
the flux intercepted by the detectors originated from in- 
elastic scattering in the target cryogen. Figure 1 shows 
charged particle spectra obtained with dedicated low- 
intensity runs and measured by drift chambers in front 
of the detectors. The dominant background was quasi- 
elastic scattering from target windows, separately mea- 
sured using an equivalent aluminum target and computed 
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to be 1.8 ± 0.2% (^He) and 0.76 ± 0.25% (^H). 

An electron must give up more than 19 MeV to break 
up the ''He nucleus and undergo quasi-elastic scatter- 
ing off nucleons. Figure 1 shows that the quasi-elastic 
threshold lies beyond the edge of the detector. A limit of 
0.15±0.15% on this background was placed by detailed 
studies of the low- intensity data. For -'H, the 7r° thresh- 
old is beyond the extent of the plot; direct background 
from inelastic scattering is thus negligible. 

Background from rescattering in the spectrometer 
apertures was studied by varying the spectrometer mo- 
mentum in dedicated runs to measure inelastic spectra 
and to obtain the detector response as a function of scat- 
tered electron energy under running conditions. From 
these two distributions, the rescattering background was 
estimated to be 0.25±0.15% (''He) and 0.10±0.05% (^H). 

For each source of background, a theoretical estimate 
for Apy was used, with relative uncertainties taken to 
be 100% or more to account for kinematic variations and 
resonance contributions. The resulting corrections and 
the associated errors are shown in Table III. Upper limits 
on rescattering contributions from exposed iron in the 
spectrometer led to an additional uncertainty of 5 ppb. 

Nonlinearity in the PMT response was limited to 1% in 
bench-tests that mimicked running conditions. The rela- 
tive nonlinearity between the PMT response and those of 
the beam intensity monitors was < 2%. A nuclear recoil 
technique using a water-cell target [11] was used to de- 
termine the scattering angle 0iab, thus keeping the scale 
error on (Q^) due to 9\ah to be < 0.2%. The acceptance 
correction K accounted for the non-linear dependence of 
the asymmetry with Q^. 

The beam polarization, Pb, was continuously moni- 
tored by a Compton polarimeter; results, averaged over 
the duration of each run, are listed in Tab. HI. Re- 
dundant cross-calibration of the recoil Compton elec- 
tron spectrum restricted the relative systematic error to 
w 1%. The results were consistent, within systematic 
uncertainties, with those obtained from recoil Compton 
photon asymmetries, and with dedicated measurements 
using M0ller scattering in the experimental hall and Mott 
scattering at low energy. Throughout the asymmetry and 
background analysis, blinding offsets were maintained 
on both results. These offsets, which were significantly 
larger than the respective statistical errors, were removed 
only after all analysis tasks were completed. After all 
corrections: 
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phys 



-6.40 ±0.23 (stat) ±0.12 (syst) ppm, 
-1.58 ±0.12 (stat) ±0.04 (syst) ppm. 



The theoretical predictions ^^^"^ ^ns with G"* = 
were estimated using the formalism in [4] and described 
in our previous publications [11, 12]. The electroweak 
radiative corrections, calculated using the MS renormal- 
ization scheme, introduced negligible uncertainties. 
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FIG. 2; 68 and 95% C.L. constraints in tlie G% - Gm plane 
from data from this apparatus ([11, 12] and this Letter). Var- 
ious theoretical predictions are plotted with published uncer- 
tainty estimates, when available. The 1-cr bands (a quadra- 
ture sum of statistical and systematic errors) and central val- 
ues (dashed lines) from the new results alone are also shown. 



Assuming a pure isoscalar 0+ 0+ transition, A^g 
is completely independent of nuclear structure and de- 
termined purely by electroweak parameters. D-state and 
isospin admixtures and meson exchange currents are neg- 
ligible at the level of the experimental fractional accu- 
racy of ^ 3% [13]. For our kinematics (£'b=2.75 GeV, 
(Q2^ = 0.077 GeV^) we obtain A^l = +6.37 ppm. 

Electromagnetic form factors from a phenomenolog- 
ical fit to the world data at low [14] were used 
to calculate ^^S' with uncertainties governed by data 
near - 0.1 GeV^ The value used for G]^ = 0.037, 
with a 10% relative uncertainty based on new data 
from the BLAST experiment [15]. For our kinematics 
{Eb^S.lS GeV, (Q2) = 0.109 GeV^) we obtain Agg = 
— 1.66 ± 0.05 ppm. This includes a contribution from 
the axial form factor G^, and associated radiative cor- 
rections [16], of -0.037 ± 0.018 ppm. 

Comparing our results to the theoretical expecta- 
tions, we extract Gj^ = 0.002 ± 0.014 ± 0.007 at = 



0.077 GeV and G% + 0.09G|^ = 0.007 ± 0.011 ± 0.004 ± 
0.005 (FF) at = 0.109 GeV^ where the uncertainties 
in the nucleon electromagnetic form factors govern the 
last error. Figure 2 displays the combined result for these 
and our previous measurements [11, 12], taken with (Q^) 
between 0.077-0.109 GeV^. The requisite small extrapo- 
lation to a common — O.lGeV^ was made assuming 



that G|; cx Q-^ and that G 
G%, = -0.005 ± 0.019 and G] 



is constant. The values 
= 0.18 ± 0.27 (correlation 
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coefficient =—0.87) arc obtained. The results are quite 
insensitive to variations in G^, as evidenced by the neg- 
ligible change induced by an alternate fit similar to that 
in [17], where is constrained by other ^pv data. 

Figure 2 also displays predictions from selected theo- 
retical models [18-23]. Those that predict little strange 
quark dynamics in the vector form factors arc favored [22, 
23]. A global fit to all low-Q^ measurements of G% and 
G%[, similar to that performed in [17], finds that other 
measurements [6, 9, 10] which had suggested non-zero 
strangeness effects are consistent, within quoted uncer- 
tainties, with our results at = 0.1 GeV^. Due to the 
improved statistical precision and lower G^ sensitivity of 
our result, adding these other measurements in a global 
fit does not alter our conclusions. 

In summary, we have reported the most precise con- 
straints on the strange form factors at ~ 0.1 GeV^. 
The results, consistent within errors with other Apy 
measurements, leave little room for observable nucleon 
strangeness dynamics at low Q^. Theoretical uncertain- 
ties, especially regarding the assumption of charge sym- 
metry [24] , preclude significant improvement to the mea- 
surements reported here. While future experiments will 
pursue the search for non-zero strangeness at higher Q^, 
it now becomes a challenge for various theoretical ap- 
proaches to reconcile these results and enhance our un- 
derstanding of nucleon structure. 
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ported by The Southeastern Universities Research Asso- 
ciation, Inc. under U.S. DOE Contract No. DE-AC05- 
84150, and by the DOE and NSF (United States), the 
INFN (Italy), and the CEA (France). 
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